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We set up a field theory model to describe the longitudinal low energy modes in high density 
matter present in white dwarf stars. At the relevant scales, ions - the nuclei of oxygen, carbon 
and helium - are treated as heavy point-like spin-0 charged particles in an effective field theory 
approach, while the electron dynamics is described by the Dirac Lagrangian at the one-loop level. 
We show that there always exists a longitudinal gapless mode in the system irrespective whether 
the ions are in a plasma, crystal, or quantum liquid state. For certain values of the parameters, 
the gapless mode can be interpreted as a zero sound mode and, for other values, as an ion acoustic 
wave; we show that the zero sound and ion acoustic wave are complementary to each other. We 
discuss possible physical consequences of these modes for properties of white dwarfs. 


INTRODUCTION AND SUMMARY 

Consider a neutral system of particles made of quan¬ 
tum degenerate electrons, and either Oxygen {Z = 8), 
or Carbon {Z = 6), or Helium {Z = 2) nuclei, at 
mass densities p ~ (10® — lO”^) g/cm®, and temperatures 
T < {a few) X 10® K. For such high densities, average 
inter-particle separations in the system are much smaller 
than the atomic scale - hence no atoms would form, even 
if the system were cold. On the other hand, the separa¬ 
tions are much larger than the nuclear scale, hence, one 
can regard the 0,C, and He nuclei as point-like, posi¬ 
tively charged, spin-0 particles; we’ll refer them as ions 
in the present work, to emphasize that one doesn’t need 
to care about their detailed nuclear structure. 

For the 0,C, and He ions the baryon number. A, 
equals twice their charge Z, and therefore the ion mass 
densities, p ^ (10® — 10^) g/cm®, translate into the elec¬ 
tron number density that can be estimated as Jq ~ 
{p/2mp) ~ (0.1 — 0.3 MeV)®, where is the protons 
mass. The corresponding Fermi momenta are pp ~ 
(0.3 — 0.9)MeV, and hence the electrons are (nearly) rel¬ 
ativistic. Their Fermi energy significantly exceeds their 
thermal energy, as well as their two-body Coulomb in¬ 
teraction energies; therefore, the electron system can be 
regarded as a quantum degenerate gas. 

As for the ions, they’re much heavier than the elec¬ 
trons, and hence their thermal de Broglie wave-lengths 
much shorter, so that at T ~ (a few) x 10® K they generi- 
cally form a classical gas (see [1], and references therein). 
However, at lower temperatures two effects need to be 
taken into account: (a) Their two body Coulomb inter¬ 
action energies start to dominate over their thermal en¬ 
ergy; (b) Their thermal de Broglie wavelengths become 
comparable to their average separation, and they enter a 
quantum regime. As a result, below T ^ {a few) x 10® K, 
the ions may form either a classical or quantum crystal 


[2], or a quantum liquid [3], depending on concrete values 
of p and Z. 

The above described matter is believed to exist in the 
interiors of white dwarf stars (WD’s). These are stars 
that finished their thermonuclear burning process, and 
are sitting in the sky to radiate away the heat stored 
in them. They can be regarded as retired stars, slowly 
evolving from being part of luminous matter to become 
baryonic dark matter. There are abundant numbers of 
such WD’s observed in our galaxy alone, and the interiors 
of the majority of them consist of Carbon, or Oxygen, or 
a mixture of the two, with mass density p ~ (10® — 10^) 
g/cm®, although higher densities are also present in many 
of them (see, e.g., [1] and references therein). A typical 
WD starts off at temperature ~ (a few) x 10® K, and 
takes from a few to 10 Gyrs to cool down to about 10® K 
or below, becoming then directly unobservable. 

The WD cooling rate is strongly influenced by thermo¬ 
dynamic properties of the state of matter in the bulk of 
these stars, and in particular by its specific heat. The lat¬ 
ter can be calculated if one knows dispersion relations of 
low energy excitations of the substance in WD’s [1]. As 
mentioned above, depending on concrete values of the 
temperature T, density p , and ion charge Z, the ions 
could end up being in a classical gas state, or may cre¬ 
ate a classical or quantum crystal, or may condense into 
a quantum liquid state. These different substances will 
have different low energy excitations, and hence, different 
specific heats and cooling rates. Knowing accurate values 
of these rates is important, e.g., for precise determination 
of the age of the Universe. 

Irrespective of the microscopic structure of the result¬ 
ing state of the ions in WD’s, the ion system may be 
treated as a uniform substance described by an appropri¬ 
ate equation of state, at length scales much greater than 
the average inter-ion separation. We’d like to study long- 
wavelength collective fluctuations in this neutral system. 
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Furthermore, we’ll be solely interested in the longitudi¬ 
nal low energy excitations for reasons that will become 
clear below. In particular, we’d like to understand the in¬ 
terplay between the zero sound mode, that typically ex¬ 
ists in degenerate fermionic systems, and the ion acoustic 
wave that is usually present in a neutral plasma. What 
we’ll show is that these two modes are complementary to 
each other: when one of these modes is present the other 
one is absent. However, one of them is always present. 
We will also show that the cores of O and C WD’s with 
p ~ (10® —10^) g/cm® support only the ion acoustic wave; 
the zero sound mode could be present in those WD’s in 
a relatively narrow spherical shell, away from the cores if 
density in those domains is p < 10® g/cm®. 

In Section 1 we will set up a prototype effective 
field theory model to calculate dispersion relations for 
the low energy excitations. Our calculations will focus 
on temperatures below (a few) x 10^ K. The two- 
body Coulomb energy of the ions is of the order of, 
(10^ — 10®) eV = (10® — 10®) K. While the Coulomb 
interactions are screened by the electrons, the screening 
length is greater than the average inter-ion separation, 
thus leading only to a small reduction, at the level of 
10% or so, of the unscreened two-body Coulomb interac¬ 
tion energy [4]. Therefore, for T < {a few) x 10^ K, the 
screened Coulomb energy dominates over the ion ther¬ 
mal energies, by at least an order of magnitude. Hence, 
we neglect the thermal effects. Likewise, finite temper¬ 
ature effects are negligible for the calculation of disper¬ 
sion relations in the fermionic part of the system, since 
TIJq <C 1. Thus, in Section 1, we formulate an ef¬ 
fective field theory at zero temperature, to calculate the 
dispersion relations for low energy modes. One can then 
use the standard formalism of finite temperature statis¬ 
tical mechanics to evaluate the effects of the dispersion 
relations on the values of the specific heat. This is a 
self-consistent procedure, as long as TjJ^' <C 1. 

In Section 2 we recover all the results of Section 1 in 
a Coulomb gauge. The advantage of the latter is that 
the properties of the longitudinal collective modes are 
captured by the phase of a scalar field describing the 
system of spin-0 charged ions. It is in this section that 
we show that the O and C WD’s will dominantly support 
the ion acoustic wave in their cores, even before the ions 
turn into a crystal state. 

Section 3 is dedicated to the He WD’s. There are 
WD’s that have a helium core due to removal of matter 
and energy from them by their binary companions (see 
discussions and references in [5]). Among these, further¬ 
more, there is a very small subclass of the dwarf stars for 
which the temperature at which the de Broglie wave¬ 
lengths of the nuclei begin to overlap, is higher than the 
would-be crystallization temperature. Then, right below 
Tc, the quantum-mechanical uncertainty in the position 
of the charged nuclei is greater that the average inter- 
nuclear separation. This is exactly opposite to the crys¬ 
tallized state where the nuclei would have well-localized 
positions with slight quantum-mechanical fuzziness due 


to their zero-point oscillations. 

It was argued in Refs. [3, 6] that such a system, instead 
of forming a crystalline lattice, would condense owing to 
the quantum-mechanical probabilistic attraction of Bose 
particles to occupy one and the same zero-momentum 
state, and leading to a quantum liquid in which the 
charged spin-0 nuclei would form a macroscopic quan¬ 
tum state with a large occupation number the charged 
condensate. 

The dispersion relations for the quasi-modes of the 
charged condensate were derived in [3, 6]. These re¬ 
sults were obtained by using the unitary gauge, and the 
Thomas-Fermi (TF) approximation for the electrons. Us¬ 
ing these results the cooling of the He core WD stars 
was studied in [5, 7], with the conclusion that they’d 
cool faster due to condensation, and that this prediction 
could be tested if a large enough sample oiHe core WD’s 
existed. 

However, it was subsequently shown in [8], the TF ap¬ 
proximation combined with the unitary gauge, misses one 
gapless mode, which could affect the cooling calculation. 
In particular, instead of using the TF approximation. 
Ref. [8], took into account the one-loop fermion effects, 
and unveiled the gapless mode. Furthermore, it was also 
shown in [8], that the gapless mode makes a negligible 
contribution to the specific heat at relevant temperatures 
and densities for He WD’s, and therefore the predictions 
for the fast cooling rates obtained in [5, 7] remain un¬ 
changed. 

In the present work we’ll confirm the results of Ref. [8] 
on the existence of the gapless mode in charged conden¬ 
sate, by doing calculations in both unitary and Coulomb 
gauges. The latter calculation turns out to be simpler, 
and does not require going beyond the TF approxima¬ 
tion to establish the existence of the gapless mode (the 
one-loop result is needed though in the Coulomb gauge to 
calculate the imaginary part of this mode). Furthermore, 
the Coulomb gauge calculations lead us to the arguments 
that the gapless mode found in [8] is an ion acoustic wave, 
that is present in charged condensate even when the finite 
temperature effects are ignored. 

A brief outlook is presented in the last section of the 
paper. 

1. A PROTOTYPE MODEL 

Consider a high density medium composed of nuclei 
and degenerate electrons, as discussed in the previous 
section. Depending on the temperature and/or compo¬ 
sition, the nuclei could be in a state of plasma, crystal, 
or charged condensate. In all the three cases Lorentz in¬ 
variance is broken by the medium; for a crystal the rota¬ 
tion and translation symmetries are also broken to their 
discrete subgroups. We consider this medium at length 
scales larger than the average inter-particle separations, 
and would like to study long wavelength longitudinal col¬ 
lective modes. We use an effective field theory technique 
to write down the Lagrangian that respects symmetries 
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of the system at hand. In this section we discuss a pro¬ 
totype Lagrangian. 

The ions are much heavier than the electrons, and at 
low momenta the effects they produce are presumed to 
be modeled by introducing in the Maxwell Lagrangian 
the “electric” and “magnetic” masses, denoted by toq 
and m^, respectively (electrons and their effects will be 
accounted for a bit later) 

= {E] - -f mlAl - mi^A]) , (1) 

where Ej and , j = 1,2,3, are the components of the 
electric and magnetic helds, and unless stated otherwise 
we use the units c = h = 1 , for which the vacuum values 
of the dielectric constant and magnetic permeability are 
set to unity. 

The presence of the “electric” and “magnetic” mass 
terms in ( 1 ) captures the phase of the spin -0 ionic mat¬ 
ter, as will be more clear below. This phase is what’s 
capturing the relevant low energy physics of longitudinal 
modes. In reality, the parameters ttiq and would be 
scale dependent quantities, and should account for renor¬ 
malization of the vacuum values of the dielectric constant 
and magnetic permeability, due to the quantum effects of 
the ions. However, we are interested in the low momen¬ 
tum and frequency limit, where these parameters can be 
approximated by constant mo and m-y. The scale de¬ 
pendent renormalization of these quantities due to the 
electrons will be explicitly included below. We’ll show 
that the outlined approximations give a reasonable de¬ 
scription for longitudinal modes in the plasma and crys¬ 
tal (Sections 1 and 2), however, they break down for the 
charged condensate (Section 3), where the explicit scale 
dependence of the “electric mass” cannot be ignored. 

In particular, the prototype model allows for a simple 
and clear description of the conditions for the existence 
of a zero sound and ion acoustic wave in the degenerate 
plasma and crystal; with a certain modification, derived 
in Section 3, this model can also describe the collective 
fluctuations of the charged condensate. 

The Lagrangian (1) is not gauge invariant. The only 
physical meaning of the latter statement is that ( 1 ) de¬ 
scribes more degrees of freedom than the Maxwell theory. 
Local gauge invariance can always be restored at the ex¬ 
pense of introducing new fields, and therefore, it is a re¬ 
dundancy of the description (although a convenient one). 
In particular, (1) can always be regarded as a gauge-fixed 
version of a gauge invariant Lagrangian obtained from (I) 
by the substitution, with the in¬ 

variance transformation, = 9 ^ 7 , 5a = 7 , where 7 
is a gauge transformation parameter, and /i = 0,1,2,3. 
The phase field a makes explicit the presence of the de- 
gree(s) of freedom beyond the two transverse states that 
can be attributed to the gauge held when a nonzero 
a is retained. 

Conversely, (I) can be regarded as the Lagrangian in 
the so-called unitary gauge, 0 = 0. We use this gauge in 
the present Section, while in Section 2 we restore back 


a nonzero a, and use instead the Coulomb gauge for the 
gauge held, Aj = 0. The results in the two gauges will 
naturally be the same, but the two derivations are differ¬ 
ent, each having its own advantages for understanding of 
the hnal results. 

Irrespective of the gauge choice, the system described 
by (I) contains an extra longitudinal degree of freedom, 
in addition to the two transverse modes of a photon. The 
longitudinal mode can be thought of as a collective low 
energy excitation of the charged ion background. Since 
we have not introduced yet the neutralizing electrons in 
( 1 ), the spectrum of excitations is gapped by the param¬ 
eter m^. As long as the latter is smaller than the inverse 
inter-particle distance, all three gapped modes can still 
be meaningfully described by the low energy Lagrangian. 

We now introduce the electrons. Instead of writing an 
effective Largangian for them, we use the fundamental 
description in terms of the Dirac theory. This is justi¬ 
fied: the ions are much heavier than the electrons, and 
at energy scales below the ion mass their collective lon¬ 
gitudinal mode is captured into an effective Lagrangian 
( 1 ), while the electrons can be kept fundamental as long 
as they are weakly interacting, and as long as their loop 
effects will in the end be taken into account (see below). 
Hence, the total low-energy prototype Lagrangian for the 
electron-ion system reads as follows: 

C = Ca+Cf^ Cf = i’iilP - me)i^ ■ ( 2 ) 

Here we have introduced the chemical potential for the 
electrons fie via the usual prescription on the covari¬ 
ant derivative Dq = Dq — ifie- We will also package 
the electric and magnetic fields in a Lorentz-invariant 
Maxwell form, with the Lorentz-breaking effects of the 
nuclei summarized by mo and m-y] this is pending addi¬ 
tional Lorentz-violating effects to arise due to the elec¬ 
trons. To account for the latter, we look at the contribu¬ 
tion of the electrons to the Lagrangian at the one-loop 
level, namely via the photon self-energy H^'^. Then, the 
Lagrangian density for the fluctuations in the quadratic 
approximation is given by 

^ = -W. + IrnlAl - im^Af + . (3) 

Note that the inclusion of a one-loop expression for H^'^ 
in the effective Lagrangian (3) to determine the disper¬ 
sion relations, yields the same results for the disper¬ 
sion relations, as if they were deduced from the poles 
of Green’s functions in which the one-loop bubble dia¬ 
grams have been resummed (such a resummation is often 
referred as the random phase approximation). 

It will be convenient to work in the momentum space, 
and use the Fourier transform of H^^ that we denote by 
. Due to the conservation of the fermion number, 
and rotation symmetry, H'^'^ can be expressed in terms 
of two functions, n(u;,A:) and n-'-(w,fc), where k = \k\, 
and takes the form 

H ^H 
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If we decompose the photon into transverse, longitudinal, 
and time-like components, ^ Aq, the Lagrangian 

for the transverse modes decouples from that of the longi¬ 
tudinal and time-like components. In momentum space, 
we have 




. 2 ^ A± 


where 


M = 


C\ = \(A„ . 


_ ( {l + +ml Lok (I -I- 

ujk{l + §) uj'^{l + ^)-r. 


( 5 ) 

( 6 ) 

( 7 ) 


The dispersion relations of the two transverse modes 
are clearly given by 

uj'^ — k'^ — (w, k) — rn^ = 0 . (8) 

The dispersion relations for the remaining modes are 
given by the zeros of the determinant of the matrix M, 


det IM = (toqo;^ — (I -f 


n(u;, k) 

*2 


mgm^ = 0. 


( 9 ) 


If we were to neglect the contribution of the electrons, 
i.e., to set n(a;, k) = 11-^ (w, k) = 0, we would have found 
from (9) one massive longitudinal mode 

0 ) = , ( 10 ) 


alongside with two massive transverse modes described 
by (8). Let us see now how the effects of the electrons, 
encoded in n(a;,fc) and n-^(w,A:), modify this spectrum. 

The expression for the non-vacuum contribution to the 
one-loop self energy of the photon due to electrons at 
finite chemical potential is standard and can be found, 
for example, in [9, 10]. In the approximation of small w 
and k as compared to the Fermi energy and momenta, 
and for u) ^ vpk, this expression is given by 




In 


OJ ■ 


k vf 


OJ — kvF 


( 11 ) 


where the Fermi velocity wp is defined as vp = kp/me- 
The usual Debye screening mass is denoted by rus 

2 2 

m2 = n(w = 0,fc^0) = (12) 

The above expression is for relativistic electrons, while in 
the non-relativistic case one has: = e^Wefep/Tr^. 

Using equation (11), we can now check for both mas¬ 
sive and massless longitudinal modes. For massive modes 
we take k —>■ 0 while keeping oj hnite. We find, 

oj'^{k —>■ 0) = -I- . (13) 

Thus the Debye screening mass contributes to the longi¬ 
tudinal mode of the massive photon. 


To check for a massless pole, we set oj = xvpk, and 
also introduce the following notations: 




.,2 ’ 


6 ^ = 


ml 


(14) 


We then take the fc —)■ 0 limit while keeping x fixed. Then 
det M = 0 corresponds to, 


- a ^) (^1 - 1 1 „ 


X — 1 


= b\ 


(15) 


The right hand side of the above expression is clearly 
both real and positive. Thus, in order for a solution to 
this expression to exist, there must be some value of x, 
either real or complex, for which the left hand side is 
both real and positive. We’ll investigate these solutions 
below. 


1.1 Zero Sound 


Consider real x that is also large, x 1. Then equa¬ 
tion (15) takes the form 


{x 


2 




= ¥ 


and the solution is 


X 


2 


5a^ — 3 
15&2 -H 5 ■ 


(16) 


O ur assu mption of a; 1 is valid as long as a 

V1 + 352 . In that case, ~ a^/(l + 36^), and the dis¬ 

persion relation for the massless mode is 


yiTsufT^ 


(17) 


Thus, the velocity of the massless mode (both phase 
and group) is given by vq/VI -f 352. niode exists 

when Vo 1$> because only in that case 

a; ^ I. Furthermore, in our system of electrons and ions, 
vp ^ Vs, and (17) can be approximated by the dispersion 
relation w ~ rig fc. 

Note that the slope of the linear dispersion relation of 
the mode (17) is larger than vp since a; » 1, and thus the 
Lj{k) line never intersects a region in which the near-the- 
fermi-surface modes live (see Fig. 1). Hence, dumping of 
the mode into electron-hole pairs is negligible. 

This excitation corresponds to the zero sound mode 
that exists in interacting degenerate Fermi systems, e.g., 
in repulsively interacting Fermi liquids and gases alike 
[11]. More will be said on this in the next section. 


1.2 Ion Acoustic Wave 

Let us now turn to the opposite regime, when x could 
be complex but both its real and imaginary parts are 
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FIG. 1. Spectrum of longitudinal collective modes and 
particle-hole excitations. 


small, i.e., x ^ 1 -I- zl. In this case, we need to take 
into account that the log entering (15) is a multivalued 
function in the complex plane of its argument. Then, it is 
instructive to start by proving that the equation (15) has 
no solutions on the first Riemann sheet of the complex 
x-plane if |x| < 1. 

To proceed with the prove, we introduce the definition 


W{x) = (x^ 



X -I-1 \ 
_x — 1_ / 


(18) 


Let us hrst consider the case that x G R and — 1 < x < 1. 
For this range of values we can use the following integral 
representation. 


W{x) 



(19) 


The epsilon prescription is given to recover the correct 
imaginary part of n(x), 

a;2 _ c + l 

ImlF(x) = - / dzzTrS{z — x), 

2 J-I 

= \x{x^-a^). ( 20 ) 

We see that for these values of x, W{x) is complex, with 
the exception of x = a. However, at x = a, we have 
lT(x) = 0 and thus these values of x do not contain a 
solution of (15). 

For all other values of x, real or complex, we use the 
following representation, 

lF(x) = (x^ - a^) f dz (21) 

Jo z — x 

For both real x and purely imaginary x, we see that W (x) 
is real. However, only when a > 1 can W (x) be positive 
and thus solve (15). 


Finally, we check for complex x. Setting x = a + i/?, 
with (T,/3 G R and a,f3 ^ 0, and using the representa¬ 
tion (21) we find through some basic algebra that either 
ImlF(x) ^ 0 or ReFF(x) < 0. Thus only in the case 
that a > 1 is W(x) both real and positive. Otherwise, 
we cannot solve (15) and no massless pole exists on the 
first Riemann sheet. This completes the proof. 

Thus we look for solutions on the second Riemann 
sheet, and rewrite the equation (15) in the following form 


- a-) (l - I (l„ ‘ 


1 X 


1 — X 


- ZTT = 6^ . (22) 


In the approximation of small x, this has a solution: 

X zz \/■ (23) 

Thus the physical mode at hand has the following dis¬ 
persion relation 

vl + k . (24) 


The sign of the imaginary parts above is such that the 
mode is on the second Riemann sheet of the complex 
energy plane. Thus, it’s a resonance, with a momentum 
dependent width. In what follows we will argue that this 
is nothing but the ionic acoustic wave. 

Consider the toq —>■ oo limit; the above mode does exist 
in this limit and its dispersion relation reads 

2 

Lo = (vs — i-^)k . (25) 

dup 

It is instructive to go back to the Lagrangian and cal¬ 
culate the exact k dependence of the mode, not only its 
low k limit. The uiq —>■ oo limit settles us in the regime 
Aq = 0. The Lagrangian then takes the form: 

C = . (26) 

This Lagrangian describes four degrees of freedom: two 
massive helicity-1, one massive helicity-0, and one mass¬ 
less helicity-0 states; the latter is the mode with the 
dispersion relation (25). To see all this we proceed by 
splitting the transverse and longitudinal parts of Aj as 
follows: 

A,=Af + ^-d^. (27) 

Upon this substitution the Lagrangian splits into two sep¬ 
arate parts, one for Af, 

\{d,Aff + \Af^Af - f + lA^n^^Af , 

and one for tt. 


1 / didjir V 

2 \ J 


1 

2 




{dtTr)Il{dtTT) 

m 2 


( 28 ) 
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It is straightforward to see that the dispersion relations 
for the two modes in Aj~ are determined by: 

uj]_ = +n-‘-(Lu±,k), (29) 

while the dispersion relations that follow from the La- 
grangian (28) are determined by the equation 




k^ + n(k,oj) 


(30) 


In the absence of 11 (i.e., if the electrons are assumed to 
be frozen) we get from (30) 


2 2 
w = 


Ze^Jo 

iriH 


(31) 


The latter is a plasmon dispersion of a charged ion gas. 
In this limit of “frozen” electrons (11 —)■ 0) there is noth¬ 
ing to compensate for a displaced change of ions upon 
their perturbation, and therefore the mode is gapped. 
Note that in the 11 = 11*^ = 0 limit the longitudinal and 
transverse modes have the same gap, w(fc = 0) = (jj±{k = 
0 ) = m~^, consistent with the rotation symmetry. 

Once the dynamical electrons are returned back, then 
the Lagrangian (28) and the dispersion relation (30) de¬ 
scribe two modes, one massive and one massless. Let us 
start with the massive mode, taking fc —)■ 0 in (30) with 
fixed and nonzero w, we find 11 ~ —m^Upfc^/(3w^), and 
get a solution = m'^+m'^Vp/3, which agrees with (13). 
The latter sets the gap for the longitudinal photon. 

For small oj and fc, on the other hand, we get a mass¬ 
less mode as follows: in this regime the real part of 11, 
Rell ~ rUg, and therefore, the real part of the corre¬ 
sponding dispersion relation reads 


rrv^k'^ 


k^ 




mt 


(32) 


This dispersion relation coincides with the toq —>■ oo limit 
of the one in (24), up to the imaginary part that can 
also be straightforwardly recovered by keeping it in 11. 
Therefore, this mode is an ion acoustic wave; sometimes 
it’s also referred as an ionic sound [12]. It corresponds 
to a collective excitation of the ions and neutralizing fast 
electrons, producing a gapless longitudinal mode. 

A few comments before we turn to the next section. 
We reiterate that the propagator in the considered regime 
has no low energy poles on the first Riemann sheet, nei¬ 
ther real nor complex ones. The obtained complex pole 
describing the ionic sound is on the second sheet. The 
function that the propagator is proportional to 


1 

( x ^ — 0^)11 — ’ 


has a branch cut and exhibits a “bump-like” behavior 
for real values of energy and momentum. Why does this 
mode have a resonance-like nature? The reason is simple: 
the phase velocity of the mode equals to xvp, that is less 


than vp- The range of allowed energies of the electron- 
hole pairs near the Fermi surface (depicted in Fig. 1) 
have the upper boundary with the slope near the origin 
that equals to vp. Thus, a wave with the dispersion w ~ 
xvpk would be damped due to near-the-Fermi-surface 
excitations, as long as a; < 1. 

This also explains why this mode does not exist for 
X ^ 1; if it existed it would have been faster than the 
near-the-Fermi-surface excitations, which then would not 
be able to catch up with the mode to give rise to a neutral 
acoustic wave; instead, in the a; ^ 1 regime only the zero- 
sound mode exists. Thus, the zero sound and ion acoustic 
wave are complementary to each other. 


2. THE PHASE EIELD AND LONGITUDINAL 
MODES 

The derivations in the previous section were performed 
in the unitary gauge, i.e., the phase of the charged scalar 
field describing collective motion of the ions was gauge- 
fixed to zero. This was reflected in the gauge non¬ 
invariant terms in Lagrangian (1). The latter was re¬ 
garded as a gauge-fixed version on a gauge invariant La¬ 
grangian, obtained from (1) by a substitution, —> 

with the unitary gauge corresponding 
to a = 0. In this section we’d like to keep nonzero a, 
but we use instead the Coulomb gauge Aj = 0. This 
will enable us to obtain and understand the results of the 
previous section in a more clear way. 


2.1 A Warmup Example 

We start with a field theory containing a complex 
scalar field, degenerate fermions, and an Abelian gauge 
field that couples to both the scalar and fermions in a 
conventional way. The charged scalar field is thought to 
model properties of charged nuclei, while fermions model 
the electrons, and the Abelian gauge boson models a pho¬ 
ton. For illustrative purposes, we consider the case when 
the U{1) gauge symmetry is spontaneously broken, and 
the gauge field acquires the mass term as a result of this 
breaking. Furthermore, we consider the parameter space 
for which the scalar field radial mode is a heavy state 
that can be decoupled from the rest of the fields. In this 
approximation, and in the Coulomb gauge, the relevant 
part of the Lagrangian reads as follows: 

+ 2^0 (—A -I- ft) Aq . (33) 

The first term is the gauge boson mass term, while the 
second term contains a part of the Maxwell Lagrangian, 
as well as the term generated due to one-loop renor¬ 
malization of the polarization operator via the fermion- 
antifermion pair. These are all the terms in the La¬ 
grangian that contain Ag and a. 
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Integrating out Aq we get for 7 ^ 0 


As before, we integrate out Aq to get 


An — 


m 


(-A + : 


^doa. 


n 


(34) 




(-A + . 


n 


-doa. 


(39) 


(we neglect an irrelevant zero mode of the inverse of the 
operator multiplying dga on the right hand side of (34)). 
Substituting (34) back into the Lagrangian (33), and cal¬ 
culating the dispersion relation for the remaining field a, 
we obtain the following expression 


-pn(a;,/c) 

-I- n(u;, k) + m? 


= . 


(35) 


Let us now study this relation in two different regimes, 
a; ^ 1 , and a; <C 1 , where as before we introduce u) = 
xvpk. For the a; » 1 case, one gets 11 ~ 

Hence, the solution to (35) reads: 


x^ = 


1 


Up -I- im?/ml 


(36) 


When ms ^ m and vp «C 1, we get the condition, a; ^ 1, 
required by the approximation made. As per the argu¬ 
ments of the previous section, this is the dispersion rela¬ 
tion of the zero sound mode 


UJ = 


_WF_ 

y^Vp + im?-/ml 


k. 


(37) 


We see that in the Coulomb gauge the zero sound mode is 
described by the phase a. Since a: ^ 1, the phase/group 
velocity of this mode is greater than vp, hence, this mode 
experiences no damping in the approximation we use. 

In the opposite regime, a; <C 1, the gapless mode does 
not exist, as shown by reductio ad absurdum: assuming 
that a; <C 1 , one gets A ~ ml{l + O{x)), and the solution 
X = Up ^ ^/(l + m?Irn?s) > 1) which contradicts the initial 
assumption that x were small. 

The lack of the ion sound wave in this model has a 
reason: we used the approximation when the Lorentz- 
invariant vacuum condensate of the charge scalar is the 
only source for the mass term (the first term in (33)); this 
implies that the number of dynamical scalars that can 
fluctuate is zero, in the approximation used (the Lorentz 
invariant Higgs vacuum has zero scalar number). Hence, 
one should not expect to have the ion acoustic wave, in 
the approximation when the number of ions that can fluc¬ 
tuate is ignored (the charge that neutralizes the electrons 
in this approximation is not dynamical). 


2.2 A More Realistic Model 

We now consider the case when the “electric” and 
“magnetic” masses of the photon are different, and the 
phase is explicitly kept, while the Coulomb gauge is as¬ 
sumed for Aj , as in the previous Subsection. The relevant 
part of the Lagrangian in this gauge then reads: 

\ + ^Aq (-A -f H) Aq . (38) 


Substituting this back into (38), and deducing the dis¬ 
persion relation for a we obtain 


2 P + U{u),k) ^ 2 

^ k"^ + n(w, k) + TOq 


Vn = 


m: 


(40) 


In a realistic system such as a plasma and solid, mo ^ 
m-y, and the ion sound speed, vq 'C 1. 

For a; » 1, we know that H ~ —ml/(3x'^), and in this 
approximation the solution to (40) reads 


x 


2 


1 

Vp/vQ + 3ml/ml' 


(41) 


Hence, when Wg > mo and vq ^ vp, we get a: ^ 1, and 
the zero sound dispersion relation 


Vp 


\/v\lvl -f Zml/ml 


k. 


(42) 


This coincides with the dispersion relations for the zero 
sound mode found in (17). 

Unlike in the previous Subsection, however, a solu¬ 
tion also exists in the opposite limit, x ^ 1 : indeed, 
assuming that x ^ 1, one gets H ~ mg(l -I- 0{x)), 
and hence the solution x = ;^a/(1 -I- mn/ml), which 

if, vo\/(l d-mp/mg) Vp , gives rise to the following 
physical mode 


+ (43) 

V 'rni 

This is nothing but the dispersion relation for the longi¬ 
tudinal ion sound wave discussed in Section 1. 

One can give heuristic arguments for the complemen¬ 
tarity of the two modes that we’ve derived above. For 
this, we consider the regime in which Vg can be neglected 
as compared with vq and vp. This is a meaningful ap¬ 
proximation for a plasma and crystal since m^ is sup¬ 
pressed by the heavy ion mass scale, while mg is de¬ 
termined by the electron mass and chemical potential. 
Then, the phase velocity of both the zero sound mode 
and the ion acoustic wave is approximately vq- When 
Vq ^ Vp, we get X ^ 1 and the ionic sound is present; 
this makes sense - the electrons are faster than the wave 
and thus they can readily follow a perturbation of the 
ions to screen its charge and create a gapless neutral 
mode, the ionic sound. In the opposite limit, vq ^ vp, 
the electrons are slower than the would be ionic sound 
wave, therefore they cannot be effectively screening ion 
perturbations, and it makes sense that the ion acoustic 
wave does not exist. Instead, the electron fluctuations 
themselves - which are now effectively screened by the 
ambient mobile ion charge distribution - form a collec¬ 
tive mode, the zero sound. 



Having the dispersions relations derived, let us apply 
them to the system of the Oxygen and Carbon ions at 
densities p ~ (10® — 10^) g/cm®. As we’ve already dis¬ 
cussed the corresponding Fermi momentum for the elec¬ 
trons is pf ^ (0.3 — 0.9) MeV. Therefore, the electrons 
are (nearly) relativistic, with vf ~ 1 , to a good accuracy. 
This implies that the value of x given in (41) can never 
be greater than the unity, since vq < 1. Therefore, we 
conclude that the zero sound cannot be supported in the 
cores of the O and C WD’s. Instead, the dominant lon¬ 
gitudinal mode in this case is the ion acoustic wave (43), 
with the sounds speed approximated by vq ^ 1. The ex¬ 
act value of the sound speed cannot be calculated in our 
effective Lagrangian approach, however, our finding con¬ 
firms the suggestion made in [ 2 ] that a longitudinal wave, 
the one that’s similar to the longitudinal acoustic wave of 
a crystal, can be used to study the cooling of the O and C 
WD’s even when the interior is in a strongly interacting 
plasma state, and the crystal is not yet formed (i.e., from 
temperatures (a few) x 10^ K, down to (a few) x 10® K, 
after which the crystal forms). 

As to the zero sound mode, its existence requires lower 
densities. While WD’s have fairly uniform density pro¬ 
files in their bulks (excluding their “atmospheres” that 
are dominantly made of H and He, with some small frac¬ 
tions of “metals”), nevertheless the interiors are not ex¬ 
actly uniform. Typically one can have variation of den¬ 
sity amounting to a factor of 5 in the ratio of the max¬ 
imal density to average density. Therefore, it might not 
be foolish to think of relatively narrow spherical shells, 
away from the cores of low density WD’s, in which den¬ 
sities might be p < 10® g/cm®. In such shells the zero 
sound mode could be supported instead of the ion acous¬ 
tic wave. It could also be interesting to look for the 
existence of the zero sound mode in the low mass brown 
dwarfs [13], where densities are p ~ 10® g/cm®. 


3. CHARGED CONDENSATE 

We now consider in detail the spectrum of the charged 
condensate. We incorporate the effect of the background 
density of electrons through their one-loop contribution 
to the photon self-energy. An analogous approach was 
used in [14, 15] to determine the electrostatic potential 
of this system at finite temperature. More recently, this 
approach was applied in [ 8 ] to argue that the dynami¬ 
cal electrons give rise to a previously unnoticed massless 
mode. Here we confirm this finding, and argue that the 
obtained massless mode is the ion acoustic wave discussed 
in the previous section. 

A classical nonzero vacuum expectation value of a field 
$ can serve as an order parameter for the condensation 
of the helium-4 nuclei, describing a state with a large 
occupation number. Fluctuations of the order parameter 
describe the collective modes of the condensate. In the 
non-relativistic approximation, the effective Lagrangian 


for the nuclei $ and electrons ip, can be written as follows 
[ 6 ]: 

C = V{X) - + fp{iU) - m,)iP , (44) 

where 

1 ID 

A = -($*Do$-(Do$)*$)-L^. (45) 

2 2 rriH 

The covariant derivative of the scalar field is given by 
= dfj, — 2ieA^. V{X) stands for a general polyno¬ 
mial function of its argument. The coefficients of this 
polynomial are dimensionful numbers that are inversely 
proportional to powers of a short-distance cutoff of the 
effective field theory. We normalize the first coefficient 
to one so that 

V{X)=X + C2X^ + ... . (46) 

The Lagrangian is invariant under global Hs(I) trans¬ 
formations, responsible for the conservation of the num¬ 
ber of scalars. Another global I7e(l) guarantees the elec¬ 
tron number conservation. Accordingly, we can introduce 
chemical potentials for both the scalars and electrons, Ps 
and pe respectively, via the usual prescription on the co¬ 
variant derivative Dq —)■ Dq — ip. 

We could also have included a quartic interaction for 
the scalars However, as long as the quartic 

coupling is small. An ^ results will not be 

affected by this term. 

It is useful to represent the scalar as $ = and 

to work first in the unitary gauge where the phase of the 
scalar is set to zero, F = 0. At a later stage, we’ll move 
instead to the Coulomb gauge where things will become 
easier. 

When the scalar chemical potential is zero, ps = 0, 
there is a solution to the equations of motion of (44) 
with a nonzero expectation value for the scalar field 

(S) = y|^’ (47) 

where Jq is the background electron density. 

Let us consider the quadratic action around this back¬ 
ground solution. We introduce perturbations of the 
scalar field as follows: 

S(a:) = y^+V^T(x). (48) 

The factor of y/mn has been introduced for convenience. 
Let us also integrate out the electrons as was done in [8]. 
Thus the Lagrangian density for the fluctuations in the 
quadratic approximation reads 

^2 = 

— ^{djT)'^ + 2mHm.y Aqt . (49) 

Here 

= (2e)^ ^ ^ C 2 , (50) 

When the coefficient C 2 is such that mo = m.y, the disper¬ 
sion relations coincide with those of a relativistic theory, 
however, in general toq m^. 
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3.1 Spectrum in the unitary gauge 


3.2 Spectrum in the Coulomb gauge 


We now consider the spectrum of these fluctuations. 
We decompose the photon into transverse, longitudinal 
and time-like components. In addition, let us also inte¬ 
grate out the scalar mode t so that we can easily compare 
our results to those of the previous section. Again, the 
transverse modes of the photon decouple entirely. Their 
dispersion relations are given by 

— n^(a;, k) — rn^ = 0 . (51) 

For the remaining modes, we have 

4= l(ylo AL).M. , (52) 

where now 

+ + uk{l + ^) \ 

(53) 


with M = The matrix M differs from the anal¬ 

ogous matrix in our prototype model only via the ap¬ 
pearance of the term. 

The dispersion relations of these modes are again given 
by the zeros of the determinant of the matrix IM. To check 
for the presence of massive modes, we use expression (11) 
in the matrix (53) and take the fc —> 0 limit, keeping w 
finite. Setting det IM = 0, we find a massive mode 

w^(fc —)• 0) = + grrigUp . (54) 

This can be though as the mass of the longitudinal com¬ 
ponent of the photon, that now carries three degrees of 
freedom. The mass receives two contributions, one from 
the charged nuclei, m^, and another one from charged 
electrons, m^u|/3. Hence, we find that the one-loop con¬ 
tribution shifts up the mass of the longitudinal mode. 
This mode was argued to be heavy to contribute to spe¬ 
cific heat of the charged condensate at relevant tempera¬ 
tures, [6]; we see that the one loop contribution makes it 
even heavier, justifying further that it can be neglected. 

To determine the possibility of a massless pole, we set 
uj = xvpk and then take the fc —> 0 limit. Then det IM = 0 
corresponds to. 


X 


2 



X + \ 

X — 1 



(55) 


Thus we see that the charged condensate corresponds 
to our prototype model of Section 1 in the particu¬ 
lar case with a = 0. As argued above, there is a 

massless resonance solution, with the dispersion relation 
2 

w ~ (vs — z|^)fc, in agreement with [8]. This is also in 
agreement with the dispersion relation obtained in Sec¬ 
tion 1, describing the ion acoustic wave in a plasma. In 
our view, this mode has the same physical origin and in¬ 
terpretation as the ion sound wave of ordinary degenerate 
plasma [12]. We further strengthen this latter point in 
Subsection 3.3 by unveiling the hydrodynamics origin of 
the ion acoustic wave in charged condensate. 


Here, we derive the results of the previous subsection 
in the Coulomb gauge. This makes the connection to the 
calculations in the prototype model presented in Section 
2 clearer, and helps to explain the origin of the ion sound 
wave in the charged condensate. 

As before, we use — d^a. The relevant part 

of the Lagrangian in the Coulomb gauge reads as follows: 

^ (moHo - (-A -f H) Aq , (56) 

where the key difference from the model of Sections 1 
and 2 is that “electric” mass in the charged condensate 
has essential dependence on the momentum: 


Too 



4M4\ 


(57) 


This suggests that for small momenta, we always end up 
in the regime of large “electric” mass for the photon; as 
discussed in Section 1, this implies that the zero sound 
mode will be absent, but the ionic acoustic wave should 
be present, in agreement with the results of Subsection 
3.1. Let us see this explicitly in the Coulomb gauge: 

We first integrate out Aq 


An = 


rrin 


-A- 


4M'‘ 

-A 


H 


dna ■ 


(58) 


Substituting this back into the Lagrangian (56), and de¬ 
ducing the dispersion relation for a we get 


UJ 


fc2-tn(w,A:) 


TTln 


P -I- n(u;. A:) -I- mg -I- 


= Vn 


TTln 


4M<i 


A:^ . (59) 


As before, we use the notation uj = xvpk. For both k and 
w approaching zero, we only get one dispersion relation, 
and that is with, x <C 1. This dispersion relation reads: 


UJ ~ —- k . 

TTls 


(60) 


This is the ionic acoustic wave, in agreement with previ¬ 
ously obtained results. 

A final comment before we move to the hydrodynam¬ 
ics considerations. In the approach adopted above the 
fermions were treated in a one loop approximation, while 
the scalars in terms of a low energy effective field theory. 
The quantum effects of the condensed scalars were cap¬ 
tured by the order parameter Lagrangian. It is instruc¬ 
tive to check that the same results are obtained if the 
scalars are also treated via the one-loop calculations, as 
was done in [14, 15]. We briefly use this method in the 
Coulomb gauge. Then, the poles of the full propagator 
are determined by 

k^+ 11^{uj,k)+Il{uj,k) = 0 , (61) 

where the polarization operator for the bosons, n®(a;, k) 
can be calculated straightforwardly via the corresponding 
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one loop diagrams [14, 15]. The part of the loops that is 
due to the existence of the condensate reads as follows: 

2 (cj — mH)^ — 


H 


{2mH + ojY 


(w + muY — — 'mjj 


- 2 ]. (62) 


Substituting this into (61) and taking the small momen¬ 
tum limit, we get the pole at 


uj ~ —- k . 

TTls 


(63) 


The latter coincides with the result already obtained 
above (60). 


3.3 Hydrodynamic considerations 

The purpose of this subsection is to demonstrate that 
the ionic sound found in the prototype model, as well as 
its counterpart emerging in charged condensate, can be 
understood in terms of standard hydrodynamics, with 
the only difference that the charged condensate hydro¬ 
dynamics equations need to retain the pressure gradient 
term even when finite temperature effects are ignored, as 
will be shown below. 

We’re looking at a degenerate plasma of electrons and 
positively charged nuclei (or ions); as already mentioned, 
depending on temperature T, the ion mass mn and 
charge Z, the system of nuclei could be in a classical 
gas state, or may create a Wigner crystal, or may be in a 
condensed quantum liquid state. In either case, at length 
scales much greater than the inter-ion separation, the ion 
system may be treated as a uniform substance described 
by an appropriate equation of state. We’d like to un¬ 
derstand the spectrum of long wavelength longitudinal 
collective fluctuations in this system. 

Let us first ignore the temperature effects and consider 
the case when the ions are in the plasma or crystal state. 
All the hydrodynamic equations presented below for this 
case are well know, but are given to emphasize the dif¬ 
ference of these states from the charged condensate, to 
which we’ll turn by the end of the section. 

The continuity equations for the electron and ion num¬ 
ber densities - denoted respectively by Ue and nn - read 
as: 


dtrie^H + dj{ne^HVj) = 0, (64) 

while the momentum equation for the ions is 

dtVj + {vkdk)vj = (65) 

Consider small localized perturbations, small over¬ 
densities 5ne,H of the electrons and ions over their back¬ 
ground values set by Jo and Jq/Z^ respectively 

ne = Jo + 5ne, nH = \jo + 5nH- ( 66 ) 


For both the plasma and crystal state at long wavelength, 
the Poisson equation for the electrostatic potential cre¬ 
ated by these perturbations reads as follows 

AAq = —ZeSnH + eSug ■ (67) 

For simplicity we consider relativistic electrons here. The 
electron over-density can be related to the local potential 
via the Thomas-Fermi approximation: 

Ep = =—eAo + pf, (68) 

and using that rie = P%f , we Hnd 

due — ifYeAo . (69) 

Substituting this expression into (67), one gets: 

(A — 7rY)AQ = —ZeSnn, (70) 

where is the Debye screening mass 

squared, due to the electrons. We now look at the system 
of three equations (64), (65) and (70), and consider their 
linearization above the background, {hg), with 

= Jo, = ^, Ap® = 0, v'f = 0. (71) 

In the linearized equations we transform to the Fourier 
moels for all perturbations, as for instance, 

6ne{x,t) = j (Ykdu}5ne{k,uj)E^‘^*~^^^^Y (72) 


As a result get the following dispersion relation from the 
linearized system of equations (64), (65) and (70): 




fc 2 


where mZ^ = 


{Ze)\JolZ) 


mn 


(73) 


This describes a gapless collective mode; at small mo¬ 
mentum the dispersion relations reads 


uj ~ — -k, 

nia 


(74) 


that is the dispersion relations of the ion acoustic wave 
(ionic sound) for plasma, while for a crystal it describes 
a longitudinal acoustic phonon. 

Things are a bit different, however, for the case of 
charged condensate. While the continuity and Poisson 
equations remain unchanged, the momentum equation 
gets modified due to the pressure term on the right hand 
side; this can be shown from the Lagrangian formula¬ 
tion of the charged condensate given in Section 3 (the 
pressure term would certainly exist in ordinary plasma 
as a finite temperature effect, but in charge condensate 
it is nonzero even when the finite temperature effects are 
ignored.) The corresponding momentum equation reads: 


dtVj + {vkdk)vj = -djbo - ( 75 ) 

The difference is due to the fact that a gradient of the 
pressure is not negligible for perturbations in the charged 
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condensate; the respective term is kept as the first term 
on the right hand side in (75); it is determined by the 
gradient of the gauge invariant potential, 69 = — 

4:M‘^/A)6nH- 

We now combine this new momentum equation (75), 
with the Poisson (70) and continuity (64) equations, and 
easily derive the modified dispersion relation for the lon¬ 
gitudinal mode in charged condensate: 


UJ 


1 




m^k'^ 

+ ml 


(76) 


The result coincides with the real part of the dispersion 
relation for the massless mode found in [ 8 ]. It also co¬ 
incides, in the low momentum approximation, with the 
dispersion relations for the ion acoustic waves found in 
Sections 1,2, and 3 of the present work. 

An important approximation made above is the 
Thomas-Fermi method. This method will in general miss 
fermion dynamics near the Fermi surface. In particular, 
as we see, it does not capture the width of the ion acous¬ 
tic wave due to its dumping by the fermionic quasiparti¬ 
cles. However, this effect was already taken into account 
by the one-loop consideration in [ 8 ], and above in the 
present work. 


4. OUTLOOK 


White dwarf stars constitute interesting physical ob¬ 
jects, and are also important for inferring key astrophysi- 
cal and cosmological data. The theory of cooling of O and 
C WD’s is well known [2], and agrees well with observa¬ 
tions [16, 17]. While at temperatures above ~ 10^ K the 
ions can be regarded as being in a classical gas state, and 
below ~ 10 ® being in a bcc crystal state, between these 
two temperatures - from ~ 10^ AT, down to ~ 10® K 
- one is dealing with a strongly interacting plasma of 
ions that is not easy to study using fundamental electro¬ 
magnetic interactions. Instead, we used an effective La- 
grangian approach that is in general well suited to study 
long-wavelength excitations, even for strongly interacting 
systems. Our finding of the absence of the zero sound 
mode, and the presence of the longitudinal ion acoustic 
wave in the interacting plasma regime, confirms a sugges¬ 


tion made in [ 2 ] that the longitudinal acoustic wave can 
be used to describe physics of WD’s in this interval of 
temperatures where neither gas not crystal descriptions 
are valid. 

As to the zero sound mode, its existence requires some¬ 
what lower densities, p < 10® g/cm®. While WD’s have 
fairly uniform density profiles in their interiors, they’re 
not exactly uniform and are described by an adiabatic 
equation of state. Thus, the density can vary from the 
core to the outskirts of the bulk by a factor of 5 or more. 
Therefore, there might exist a subclass of low density 
WD’s, with relatively narrow spherical shells away from 
the cores, in which densities might be p < 10 ® g/cm®; 
if so then the zero sound mode could be supported in 
those domains instead of the ion acoustic wave. It could 
also be interesting to see if the zero sound mode could be 
supported in the low mass brown dwarfs stars studied in 
[13], where densities are p ~ 10® g/cm®. 

Finally, the interaction between the electrons due to 
the exchange of the zero sound mode (or the ion acous¬ 
tic wave) would be strong if the momentum transfer in 
the two-by-two electron scattering amplitude is near the 
pole of the zero sound (or the ion acoustic wave), i.e., 
is at a; ~ xvpk, with a; ^ 1 (or with x ^ 1 for the 
ion acoustic wave). This interaction, if attractive for a 
certain domain of momenta, might lead to formation of 
bound states, or loosely bound states such as Cooper 
pairs. While naive kinematical arguments suggest that 
the momentum transfer is not close to the zero sound 
pole as long as x ^ 1 (or with x <C 1 for the ion acoustic 
wave), the issue needs careful study, presumably via the 
Schwinger - Dyson equation, to see if such interactions 
can be used to trigger Cooper instability and produce a 
gap of a physically meaningful magnitude. 
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